Analogues of D-branes in Bose-Einstein Condensates 
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We demonstrate theoretically that analogues of D-branes in string theory can be realized in 
rotating, phase-separated, two-component Bose-Einstein condensates and that they are observable 
using current experimental techniques. This study raises the possibility of simulating D-branes in 
the laboratory. 
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I. INTRODUCTION 

String theory is the most promising candidate for pro- 
ducing a unified theory of the four fundamental forces 
of nature. Dirichlet (D-) branes, non-perturbative soli- 
tonic states of string theory, have been the most funda- 
mental tool for studying non-perturbative dynamics in 
string theory. They are characterized as hypersurfaces 
on which open fundamental strings can terminate with 
the Dirichlet boundary condition |l|. D-branes are dy- 
namic objects and their collective motion is described 
by the Dirac-Born-Infeld (DBI) action in the low-energy 
regime 0, which is a nonlinear action of the scalar field 
(corresponding to the transverse position of the D-brane) 
and the U(l) gauge field 0,13 • Since the discovery of D- 
branes, string theory has developed in conjunction with 
the study of D-branes [|| . 

Several years after the discovery of D-branes, a D- 
brane-like soliton was found in field theories such as the 
nonlinear sigma model (NLcrM) and gauge theory Q , 
in which vortex strings terminate on a domain wall be- 
tween different vacua. Here, the domain wall can be iden- 
tified with a D-brane in a sense that its collective motion 
is also described by the same DBI action as that in string 
theory, while vortex lines attached to it are analogous to 
fundamental strings because their endpoints are electri- 
cally charged and identical to solitons known as "Blons" 
in the DBI action [H, H|. From this property, the au- 
thors of Ref. @ called it "D-brane soliton" . These the- 
ories thus offer simplified models for studying D-branc 
dynamics that are much easier to analyze than in full 
string theory. All possible solutions of the wall-string 
composite soliton have been classified and constructed in 
more general sigma models and gauge theories [HI [H| • 

Bose-Einstein condensates (BECs) of ultracold atoms 
are extremely flexible systems for studying solitons 
(or topological defects) since optical techniques can be 
used to control and directly visualize the condensate 
wave functions [HJ. Interest in various topological de- 
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fects in BECs with multicomponent order parameters 
has been increasing; the structure, stability, and cre- 
ation/detection schemes for monopoles I13l-[l5ll . three- 
dimensional Skyrmions fl6l - fl8j , and knots |19| have been 
discussed. In the present study, we consider a three- 
dimensional composite soliton in two-component BECs 
that consist of a domain wall and quantized vortices ter- 
minating on the wall, as sketched in Fig. 1(a). Spe- 
cific examples of the system include a BEC mixture of 
two-species atoms such as 87 Rb- 41 K [H or 85 Rb- 87 Rb 
[2lj . where the miscibility and immiscibility can be con- 
trolled by tuning the atom-atom interaction via Fesh- 
bach resonances. Here, the domain wall is referred to as 
an interface boundary of phase-separated two-component 
BECs. Although the interface has a finite thickness, the 
wall is well-defined as the plane in which both compo- 
nents have the same amplitude. Since a description of 
two-component BECs can be mapped to the NLcrM by 
introducing a pseudospin representation of the order pa- 
rameter |22j | , the resultant wall-vortex composite soliton 
corresponds to the D-brane soliton described in Ref. [6j. 
Here, the domain wall of the two components can be iden- 
tified as a D-brane soliton because it has a localized U(l) 
Nambu-Goldstone mode which can be rewritten as U(l) 
gauge field on the wall 6], which is a necessary degree 
of freedom for the DBI action of a D-brane. In addition, 
vortex lines attached to the domain wall are identified 
as fundamental strings since their endpoints are electri- 
cally charged and identical to solitons called Blons in the 
DBI action. The primary differences from the D-branc 
soliton in Ref. |6| are only that our system is nonrela- 
tivistic and does not have supersymmetry. We find that 
these composite solitons are energetically stable in rotat- 
ing, trapped BECs and are experimentally feasible with 
realistic parameters. 



II. SYSTEM 

The order parameter of two-component BECs is = 
(4 , i,\E , 2 ), where ^ = y^nje 103 (j = 1,2) are the macro- 
scopically occupied spatial wave function of the two com- 
ponents with the density rij and the phase 9j . The order 
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FIG. 1: D-brane soliton to which strings (vortices) attach in two-component BECs. (a) Schematic illustration of the wall-vortex 
soliton configuration viewed on a length scale larger than the domain-wall width and the vortex-core size. The two-component 
BECs (z > 0) and \l/2 (z < 0) are separated by the domain wall in the z — plane. The boundary condition is given by 
S z — > 1 ( — 1) for z — > oo (— oo). We assume that vortex lines are straight and perpendicular to the wall, (b) The simplest 
D-brane soliton in which a single vortex along the z-axis for z < is connected to the domain wall. The spin texture (denoted 
by arrows) indicates that the connecting point can be identified as a monopole, as seen in the enlarged view of the region near 
the connecting point, (c) The spin texture S for the solution of Eq. (|16[) with u v — 1/rj, (j>o = 0, Zo = 0, M — 1, and the 
+ sign, in the z — and y — planes. The color indicates the magnitude of S z . S = z along the vortex core (x — y = 0). 
(d) Equilibrium solutions for the phase-separated 87 Rb- 85 Rb BECs under rotation, obtained by numerical minimization of the 
energy functional for SI = 0.38uj x . The figure shows a normalized isosurface of the density difference |ni — ri2\b 3 /N — 1.5 x 1CP 3 
with the harmonic oscillator length b = ^/TiLo x /m and the total particle number N. The vortex core appears as the lighter 
region. 



parameter can be represented by the pseudospin 

S = (S x , S y , S z ) = (sin cos 0, sin^sin^, cos 6*) (1) 

with polar angle 9 = cos _1 [(rti — 7i2)/nr] and azimuthal 
angle <fi = 9 2 — 9\ as 



cos -e 
2 



, sin — e 
2 
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where nx = n\ + ni and = 9\ + 9 2 represent the lo- 
cal density and phase, respectively 22]. The simplest 
wall-vortex configuration is schematically depicted in 
Fig. Hlb), where the wall is characterized by the S z = 
plane and a vortex line in the ^ component along the 
z-axis attaches to it. Let us consider a surface in z > 
bound by the wall that encloses the end point of the vor- 
tex, shown by the dotted-curves. Then, the surface is 



mapped to a hemisphere in spin space, where the spin 
texture (S x ,S y ) = (cos <f>, sin (j>) on the wall winds once 
along the boundary. The spin below the wall varies in a 
similar manner, except along the vortex line. Thus, the 
spin texture around the end point forms a monopole-like 
configuration. 

The solution of the solitonic structure in two- 
component BECs are given by the extreme of the Gross- 
Pitaevski (GP) energy functional 
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Here, rrij is the mass of the jth component and fij is its 
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chemical potential. The system is supposed to rotate at 
the rotation frequency f2 = f2z; thus, the harmonic trap 
potential Vj — \mj{uj 2 x 2 + u 2 y 2 + uj 2 z 2 ) is modified 

by the centrifugal term as Vj = Vj — rrijVL 2 (x 2 + y 2 )/2. 
The coefficients gn, (722, and gi2 represent the atom- 
atom interactions. They are expressed in terms of the 
s-wave scattering lengths an and 022 between atoms in 
the same component and a±2 between atoms in different 
component as 



g 3 k 



TO 7 fc 



(4) 



with m~ k — mj + . The GP model is given by 
the mean-field approximation for the many-body wave 
function and provides quantitatively good description of 
the static and dynamic properties of the dilute-gas BECs 
0. 



III. MAPPING TO THE NONLINEAR SIGMA 
MODEL 

To derive the generalized NLcrM for two-component 
BECs from the GP energy functional ([3]), we assume 
?7ii = mi = m and Vi = Vi = V. By substituting the 
pseudospin representation Eq.@ of \&, we obtain [22fl 
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Vtit 



(v off - n x r) 2 + c + Cl S z + c 2 S 2 j, (5) 
where we have introduced the effective superflow velocity 

h 



V e ff : 

and the coefficients 
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co = -g-[n T (5xi + .922 + 2.912) - 4(/ji + fj, 2 )}, 
ci = ^-[n T ( gil - g 22 ) - 2(/Xi - /j 2 )], 

9 

ft 

C2 = -Z-(!9ll + 922 - 2512). 



(6) 

(7) 
(8) 
(9) 



The coefficient c\ can be interpreted as a longitudinal 
magnetic field that aligns the spin along the z-axis; it 
was assumed to be zero in this study. The term with the 
coefficient C2 determines the spin-spin interaction asso- 
ciated with S z ] it is antiferromagnetic for C2 > and 
ferromagnetic for C2 < 22j. Phase separation occurs 
for C2 < 0, which we have focused on. Further simplifi- 
cation can be achieved by assuming that Vj = and the 
total density is uniform through the relation iit = fi/g 
where g — gn = g 2 2 and fx = /ii = and that the 
kinetic energy associated with the superflow v e ff — £1 X r 
is negligible; the effects of these terms are discussed in 



the text. By using the healing length £ = %/ \/2ragn^ as 
the length scale, the total energy can reduce to 



E 



griTS, 



dr- 



J2(VS a ) 2 +M 2 (l-S 2 ) 



,(10) 



M 2 = 4|C2 



2 ,(H) 



where M is the effective mass for 5^ . This is a well-known 
massive NLcrM for effective description of a Heisenberg 
ferromagnet with spin-orbit coupling. 

The D-brane soliton by Gauntlett et al. [6J can be 
reproduced as follows. Introducing a stereographic coor- 
dinate u — (S x — iS y )/(l — S z ), we can rewrite Eq. (|T0|) 
as 



E = dr 



(1 



|2\2 



(12) 



For a fixed topological sector, vortices (a domain wall) 
parallel (perpendicular) to the z-axis, the total energy is 
bounded from below (Bogomornyi-Prasad-Sommerfield 
bound) as E > |T W | + |T V | by the topological charges that 
characterize the wall and the vortex: 



T w = M \ dr 



T v = i / dr 



u*d z u + ud z u* 

(i + M 2 ) 2 ' 

d x u*d y u — d v u*d x u 

2\2 " 



(i + M 2 ) 



Then, the equations 



d z u =F Mu = 0, (d x =F id y )u = 



(13) 
(14) 

(15) 



are satisfied, giving the analytic form of the wall-vortex 
composite solitons: u(z,r/ — x + iy) = u w (z)u v (rj), where 



u w (z) 



o ^M(z-z )-i4>a 



u v (r)) 



(i) 



\tl\iv 



(2)x 



.(16) 



The function u w represents the domain wall with wall 
position zq and phase 4>o given by (S X ,S V ); this phase 
4>o yields the Nambu-Goldstone mode localized on the 
wall. The function u v gives the vortex configuration, 
being written by arbitrary analytic functions of 77; the 
numerator represents N Vl vortices in one domain 
component) and the denominator represents N V2 vor- 
tices in the other domain (^2 component). The posi- 



tions of the vortices are denoted by r/j and 77^ 



(2) 



In 



Ref. [6j, the solution is denoted by the coordinates 
(X, if) — (tanhlog |«|, argu), where X = tanh[=FM(z — 
zo)+log |m v (?7)|]- The total energy does not depend on the 
form of the solution, but only on the topological charges 
as T w — ±M or (per unit area), and T v = 27r7V v (per 
unit length), where N v is the number of vortices passing 
through a certain z = const plane. 

Figure []Jc) shows the texture of S with the simplest 
wall-vortex configuration [Fig. [ljb)], corresponding to 
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the solution u(z,rj). A vortex exists in z < and forms 
a texture known as the lump in field theory 12311 or the 
Anderson-Toulouse vortex in superfluid 3 He [24j [. where 
the spin points up at the center and rotates continuously 
from up to down as it moves radially outward. The vor- 
tex ending attaches to the wall, causing it to bend loga- 
rithmically as z = \ogr/M [Fig. [IJc) bottom]. We can 
construct solutions in which an arbitrary number of vor- 
tices are connected to the domain wall by multiplying by 
the additional factors rj — r}j [see Eq. (fT5j)]; Fig. [2Ja) 
shows a solution in which both components have one vor- 
tex connected to the wall. In the NLcrM, the energy is 
independent of the vortex positions T/y on the domain 
wall; in other words, there is no static interaction be- 
tween vortices. 

Although analogues of "branes" have already been 
studied in the AB phase boundary of superfluid 3 He [25j . 
their exact correspondence to those in string theory has 
not been clarified. In our case, the domain wall has a 
localized U(l) Nambu-Goldstone mode and it can be 
rewritten as the U(l) gauge field on the wall, which is 
a necessary degree of freedom for the DBI action of a 
D-brane. Gauntlett et al. have shown that Eq. p^|) 
reproduces the "Blon" solutions of the DBI action for 
D-branes in string theory Q, as can be demonstrated 
by constructing an effective theory of the domain wall 
world volume with collective coordinates ZQ(x,y) and 
<fio(x,y) in u w (z). On the domain wall, Eq. be- 
comes MX = J2j=l !og(?? - r)f>) - J2j=l ~ 11 j) 
with X = zq + itfiv/M [28]. For the example shown in 
Fig. [Ua), as we travel once around infinity rj — ¥ r/e 27 ™, 
the phase angle on the domain wall world volume shifts 
as 0o — > 4>a + 27T. When we introduce the U(l) gauge 
field Aj by taking a dual as 

di<p = eijkdjAk, (17) 

the endpoints of the vortex strings are electrically 
charged particles [8j, |9j. Therefore, our domain wall can 
be identified as a D-brane on which fundamental strings 
terminate. 



IV. NUMERICAL SIMULATIONS 

To see whether the composite solitons in Eq. (| arc 
stable under actual experimental conditions, we study 
the wall-vortex soliton in a realistic setup of trapped 
two-component BECs by numerically minimizing the GP 
energy functional Eq. ([3]) (equivalently, the generalized 
NLcrM Eq. ([5])) in the three-dimensional system via 
imaginary time propagation. Then, there are additional 
contributions to the massive NLtrM: a trapping potential, 
a gradient of ut, and a kinetic energy of superflow Eq. 
dU given by the gradients of 9 and S. According to Papp 
et al. H3, we set *i (* 2 ) to a 87 Rb BEC ( 85 Rb BEC), 
and set the particle number to N = 5 x 10 4 for both 
components, the intraspecies s-wave scattering lengths 



to ai — <22 = lOOcto, and the interspecies s-wave scatter- 
ing length to ai2 — 213ao, where ao is the Bohr radius 
and where the condition for phase separation a\ 2 > 
is satisfied 26] . The experiment revealed that the value 
of 0,2 was tuned in a wide range 50cto — 900ao via Fes- 
hbach resonances [2lj . We prepare a cigar-shaped har- 
monic trap with frequencies lo x = ui y = 20 x 2n Hz and 
lu z = 5 x 2n. Rotation ft = £lz is applied to stabilize 
vortices in the condensates. 

Figure QJd) shows the isosurface of the density dif- 
ference \m — n-zl oc \S Z \ of the stationary solution for 
fl = 0.38w, representing the wall-vortex soliton corre- 
sponding to Fig. QJb) f° r trapped BECs; the regions 
of the domain wall (S z ~ 0) and the vortex core are 
clearly visible in this figure. The vortex in vf^ forms a 
coreless vortex, where its core is filled by the density of 
and transforms into a singular vortex with increasing 
distance from the domain wall. This configuration is en- 
ergetically stable since it is obtained by imaginary time 
propagation. The spin texture of this solution is almost 
identical to that in Fig. HJc), despite there being extra 
contributions in the generalized NLcrM. Since the diver- 
gent kinetic energy generated by the vortex significantly 
depletes tit to form a singular vortex core, it slightly 
modifies the potential of the spin field. 

However, when the system contains multiple vortices, 
the above effects of the extra terms become more im- 
portant. Figure [U[b) shows the equilibrium solution in 
which both components have a single vortex. The end 
point of the vortices in each component is spontaneously 
displaced from the center, while the energy is indepen- 
dent of r)j for Eq. (fT6f given by the NLcrM. This is due 
to the effective repulsion between the endpoints of each 
vortex, originated from the two energetic constraints: a 
broad distribution of vorticity near the domain wall to 
reduce the associated kinetic energy and a smooth distri- 
bution of tit to reduce its gradient. When the rotation 
is further increased, multiple vortices are generated. For 
fl = 0.80w in our parameter setting, the domain wall 
tilts to be in parallel with the rotation axis, where some 
of the vortex lines are absorbed by the wall to form a 
"vortex sheet" [27| [Fig. H^c)]. To keep the domain wall 
perpendicular to the rotation axis for the fast rotation, 
one must increase the interspecies scattering length 012, 
which decreases the interface area to reduce its energetic 
cost. Then, as shown in Fig. [UJd), the vortex endings are 
also shifted relative to each other to form an interlaced 
rectangular lattice on the domain wall. 



V. CONCLUSION AND DISCUSSION 

We have shown that an analogue of a D-brane can 
be realized as an energetically stable solitonic object in 
phase-separated, rotating, two-component BECs. This 
suggests that atomic BECs are the most promising can- 
didates for demonstrating D-brane physics in the labora- 
tory; the energetically stable D-brane solitons warrant 
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FIG. 2: D-brane to which several strings attach in two-component BECs. (a) Left: Schematic illustration of the configuration 
in which each component has a single vortex connected to the wall. The connecting points can be seen to be a monopole 
and an antimonopole. Right: The spin texture of this configuration in the z = and y = planes, using Eq. (|16p with 
u v — (rj — xo)/(r) + xo) and xo = ±2. The wall becomes asymptotically flat due to the balance between the tensions of the 
attached vortices, (b)-(d) Equilibrium solutions obtained by numerical minimization of the Gross-Pitaevskii energy functional 
for 87 Rb- 85 Rb BECs. The expression is the same as that in Fig. QJd). The rotation frequency is (b) = 0A0lj x , and (c) (d) 
fi = 0Wu x . Only for (d), oia = 473a . 



studying various dynamic phenomena, e.g., oscillation 
modes of strings and branes and nonlinear dynamics such 
as brane-antibrane annihilation, which was proposed as a 
possible explanation for the inflationary universe in string 
theory. Here, we summarize the outlook gained by this 
realization. 



A. D-brane — anti-D-brane annihilation 

Although brane-antibrane annihilation was demon- 
strated to show the topological defect creation in su- 
perfiuid 3 He (25|, a physical explanation of the creation 
mechanism of defects still remains unclear. We note 
that the intriguing experiment that mimicked the brane- 
antibrane annihilation was performed by Anderson et al. 
[3fj| . They created the configuration shown in Fig [31 
where the nodal plane of a dark soliton in one component 
was filled with the other component. By selectively re- 
moving the filling component with a resonant laser beam, 
they made a planer dark-soliton in a single-component 
BEC. It is known that the planer dark soliton in 3D sys- 
tem is dynamically unstable for its transverse deforma- 
tion (known as snake instability) [30(, which results in 
the decay of the dark soliton into vortex rings. 

In our context, this experiment demonstrated the 
brane-antibrane collision and subsequent creation of cos- 
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FIG. 3: Schematic illustration of the D-brane-anti-D-brane 
collision in the BEC system. 



mic strings, where the snake instability may correspond 
to "tachyon condensation" in string theory [31]. The 
procedure that removes the filling component can de- 
crease the distance R between two domain walls and 
cause their collision [see FigJ5]. The tachyon condensa- 
tion can leave lower dimensional topological defects after 
the annihilation of D-brane and anti-D-brane. In our 
case of the phase-separated two-component BECs, the 
annihilation of the 2-dimensional defects (domain walls) 
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leaves 1-dimensional defects (quantized vortices), which 
we would like to identify as closed fundamental strings 
because it can end on a domain wall. In string theory, 
creation of lower dimensional D-branes after D-brane an- 
nihilation was studied very well [311 ] . while creation of 
closed fundamental strings is in general difficult to deal 
with. Contrary to this, it is easy in principle in our case 
to study closed string creation in connection with tachy- 
onic fluctuations, which will be one of merits of our sys- 
tem. 

A D-brane soliton similar to ours can be identified with 
a D-brane in string theory [32lj in a tachyon effective field 
theory (known as the Minahan-Zwiebach model (33|) on a 
non-BPS D-brane. Furthermore, a brane-anti-brane an- 
nihilation was studied in the Minahan-Zwiebach model in 
[34| . There, a string connecting between the two branes 
and tachyon fluctuations were studied. In the case of our 
D-brane solitons too, one can construct a string between 
a pair of brane and anti-brane. Investigation of tachyonic 
fluctuations with or without a string between two branes 
will be reported elsewhere. 

B. Supersymmetry 

Our system does not have supersymmetry which is a 
basic ingredient in string theory. Nevertheless, our D- 



brane soliton is stable because it is topological, in con- 
trast to D-branes in string theory which are stable due 
to supersymmetry. As this concern we have two com- 
ments. One is that our system can be made supersym- 
metric without changing the bosonic part by introduc- 
ing additional fermions to the system. For example, the 
non-relativistic superstring can be realized by trapping 
the fermionic atoms in the core of vortices in a BEC 
[35j . Several studies proposed the possible simulation of 
(non)relativistic supersymmetry models using a mixture 
of ultracold fermions and diatomic bosons in optical lat- 
tices (36|. Then we expect that our D-brane soliton can 
become a BPS object preserving a fraction of supersym- 
metry, as in D-branes in string theory. In fact at least 
in the sigma model limit, our solitons reduces to a BPS 
soliton of supersymmetric theories @, 0, H3, EH ■ There 
is merit in studying further by bringing supersymmetry 
in our system. 
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